Charge-spin duality in non-equilibrium transport of helical liquids 
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Non-equilibrium transport properties of charge and spin sector of two edges of a quantum spin 
Hall insulator are investigated theoretically in a four-terminal configuration. A simple duality rela- 
tion between charge and spin sector is found for two helical Tomonaga Luttinger liquids (hTTLs) 
connected to non-interacting electron reservoirs. If the hTLLs on opposite edges are coupled locally 
or non-locally, the mixing between them yields interesting physics where spin information can be 
easily detected by a charge measurement and vice versa. Particularly, we show how a pure spin 
density in the absence of charge current can be generated in a setup that contains two hTLL and 
one spinful Tomonaga Luttinger liquid in between. 
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Introduction.- The discovery of topological insulators 
(TIs) in both two spatial dimensions (2D) and three 
spatial dimensions (3D) has recently attracted a lot of 
intereslji^— . Unlike in normal insulators, there is a gap- 
less mode appearing within the bulk gap at the edge 
of TIs which originates from strong spin-orbit coupling 
and is protected by time reversal. In 2D, a TI is also 
called quantum spin Hall (QSH) insulator since its edge 
states are one dimensional (ID) counter-propagating 
modes with opposite spin. These ID systems have been 
coined helical liquids or helical Tomonaga Luttinger liq- 
uids (hTLLs). Transport properties of hTLLs have been 
predicted and observed at the edge of HgTe quantum 
wells^i^ and proposed to also exist in InAs/GaSb quan- 
tum wells^ as well as Bi2Se3 or Bi2Te3 thin films^'^. An 
important feature of the hTLL is that spin and momen- 
tum are locked to each other. Remarkably, one hTLL has 
only half the degrees of freedom of a spinful Tomonaga 
Luttinger liquid (sTLL). Thus, two hTLLs, which natu- 
rally exist at two opposite edges of a QSH insulator, can 
recover the degrees of freedom of a sTLL. It is well known 
and has even been experimentally confirmed that there 
is spin-charge separation for a ID sTLL^. 

Therefore it is natural to ask the question how spin 
and charge sector behave for two uncoupled as well as 
two coupled hTLLs. In this Letter, we investigate the 
non-equilibrium transport properties of two hTLLs in 
a four-terminal configuration. Most interestingly, we 
find a duality relation between charge and spin sec- 
tor of two hTLLs taking into account the coupling to 
non-interacting electron reservoirs. As a physical conse- 
quence, there is a simple relation between charge current 
and spin polarization in the dual voltage configurations 
of two hTLLs (see below) . Importantly, the coupling be- 
tween two edges will destroy the simple duality relation. 
However, we can still manipulate the charge and spin sec- 
tor separately only by electric means. To demonstrate 
this, we study different scattering mechanisms between 
the two hTLLs within the non-equilibrium Keldysh for- 
malism and bosonization. Different bias dependencies 



are found for different scattering mechanisms which can 
be used to distinguish and identify them in experiments. 

Model and spin-charge duality.- We consider a QSH 
insulator in a four-terminal configuration as shown in 
Fig. [T] (a). The two edges are denoted by a = +(— ) for 
the upper (lower) edge. On each edge a, there are two 
terminals with chemical potentials ^i^a (j- = 1,2 means 
left and right lead, respectively). Two terminals on the 
same edge are connected by a hTLL of finite length L. 
The hTLL states are described by field operators V'aa 
where (a, cr) = {R, f) or (L, |) for the upper edge and 
(a, a) = {R, I) or (L, f) for the lower edge. In the middle 
region of the sample, the hTLLs at the two edges can 
mix and different types of coupling mechanisms will be 
discussed below. 

Interestingly, there are two possibilities for choosing 
the basis states of the system: the helical edge basis and 
the spin-charge basis. For the helical edge basis, the non- 
chiral boson field is defined separately for each edge with 
</?+(_) = <^i?t(;)+</'L;(T) at upper (lower) edge and the 
corresponding dual field 6'+(_) = (I>r^{i) - <l>LUt)- Here, 
4'acr {a = R, L and a =t, j,) is the standard boson field 
operator in bosonization^. This basis is suitable to study 
the current at different terminals. However, when we 
are interested in spin properties, it is more convenient to 
introduce the spin-charge basis, which is related to the 
helical edge basis by 



(1) 



where c and s represent charge and spin sector, respec- 
tively. The Hamiltonian can be written as 



H — Hq + Hy + Ht, 



(2) 



where Hq describes the hTLLs at two edges, Hy the 
coupling between the helical liquid and the leads, and 
Ht the scattering region. 
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FIG. 1: (Color online) (a) Schematic of the four-terminal 
setup. At each edge, there is a conducting channel of a hTLL 
(blue full lines correspond to spin up; red dashed lines to spin 
down). The two hTLLs are mixed in the junction region and 
different types of junctions are analyzed: (b) the short junc- 
tion with two possible single particle scattering terms; (i) spin 
conserved scattering tc and (ii) spin-flip scattering ta and (c) 
the long junction modeled by a sTLL. 



The hTLL coupled to electron reservoirs can be mod- 
eled by the so-called 5(x)-niodel"^-"'^'^ given by 



^0 = 



^0 = 
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hvf 



dx 
dx 
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g'^{x) 



(3) 



in the spin-charge basis. Here, w/ is the Fermi veloc- 
ity and g{x) is the hTLL interaction parameter {g{x) = 
go < 1 for repulsive interactions within the helical edge 
located at \x\ < L/2, and g{x) = 1 for the non-interacting 
fermions in the leads with |x| > L/2)^. The chemical po- 
tentials in the leads are naturally taken into account with 
the Hamiltonian Hy — J [dx^J.cix)'Pc + dxlJ-six)Os] 
where Oxf^ds) = -Mi,c(s)'5(a; + L/2) + fi2,c{s)S{x - L/2) 
with ^i^c = + A^i - and ^i^s = fJ-i,+ - fJ-i,^ {i = 1, 2). 
Remarkably, /ic couples to Lpc while /i^ couples to 9s- 
Therefore, the electric voltage can couple to both charge 
and spin sector. This provides us an easy way to con- 
trol charge and spin separately - in contrast to the usual 
sTLL where the electric chemical potential only couples 
to the charge sector. Moreover, we discover that there is 
a duality relation between charge and spin sector, namely 



Oc ^ (fis 



(4) 



For 



(A = c, s), the above duality relation is directly 
related to the constraint gc — I /gs discovered in Ref. [IJ] 
before. Here, we show that this relation remains valid 
even if the system is coupled to biased electron reser- 
voirs. Thus, the duality relation should be observable in 
transport properties of the system. 



What is the physical consequence of this duality re- 
lation? To answer this question, we investigate the 
total charge current and spin density of the system. 

The charge current is given by jc — e^^J^dtfc, which 

is the sum of the currents along the two edges jc = 
j+ + j_, with ja{x) — -^dtfa- The spin density can 

be defined as Ps = iJ^dx^Ps- Combining Eq. ([T]) and 

dt^a = —VfdxOa, it is evident that the spin density 
can be directly related to the charge current along the 
two edges by ps = — -i-(j+ — j-)- In the absence of 
mixing between two edges, our setup describes trans- 
port through two independent ID channels. Then, it fol- 



lows directly from previous work- --^ that (jo.q 
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where eVc 

{jo.c 



^Vc with Vc 
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the total spin density is given by (/5o. 

V+~V- _ 1 / 



Thus, the total charge current is 
}^{pi.c - P-2.c), while 

2e 



T^Vs with 

hvf ^ 



2 ~ 2eV/-"r,s A'2,s)- Importantly, it is spin 
density and not spin current that is dual to the charge 
current, which is a direct consequence of the duality rela- 
tion Q . Physically, Vc and Vs can be easily generated by 
two different voltage configurations of the four terminal 
setup, as shown in Fig. [2] In these two symmetrical bias 
configurations, we find either charge current or spin den- 
sity but no spin current. This is different for unsymmet- 
rical bias configurations where charge current and spin 
density are usually accompanied by spin current as well. 




FIG. 2: (Color online) (a). (I) and (III) show band dis- 
persions and chemical potentials of the upper edge and the 
lower edge, respectively. (II) illustrates the voltage config- 
uration p\^+ = = eVb/2 and p2,+ = P2,- ~ — eVb/2, 
yielding Vc = Vo and Va = 0; (b). Similar to (a) with 
a different voltage configuration = p2,- = eVo/2 and 
P2,+ ~ Pi.- = — eVo/2, giving Vc = and Va = Vo- Note that 
in (a) there is a finite charge current but no spin density while 
in (b) there is a finite spin density but no charge current. 

Up to now, we have discussed transport properties of 
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two hTLLs and shown a simple relation between charge 
current and spin density in two (dual) voltage configura- 
tions. In the following, we would like to take into account 
a junction structure introducing scattering between two 
hTLLs within a region of finite length d. We consider 
two different scenarios: (i) the short junction (SJ) with 
d <^ Xf and (ii) the long junction (LJ) with d ^ A/, 
where A/ is the electron Fermi wave length. For SJ, we 
can neglect the length of the scattering region and model 
it as a quantum point contact, while for LJ, we can regard 
the scattering region as a sTLL with finite length. 

Short junction case.- In the following, we concen- 
trate on the experimentally relevant regime l/-\/3 < 
go < 1—. Then, all possible one-particle and two- 
particle scattering terms will be irrelevanl^. Hence, 
we can safely treat the scattering Hamiltonian Ht 
as a perturbation. For the SJ case, two types 
of one particle scattering termsi^, which preserve 
time reversal symmetry, are taken into account (see 
Fig.^b)), given by 4 = ii,E.=± r?^"??^" sin(v^S,) + 
its X]a=± Tj^'^Tj"'^ sm{-\/TT'^ a) whcrc tj, term preserves spin 
and tg term flips spin. Here, S^- = ^/2{lpc + cfips), 
^'a = V2{ips + aOs), and 77°°' is the so-called Klein fac- 
tor. We now perform the perturbative calculation of our 
four terminal system within the non-equilibrium Keldysh 
formalis m"""— . All physical quantities can be related 
to expectation values of boson fields. Treating Ht as 
a perturbation, we can expand any physical quantity O 
{O = jcPs) in powers of e.g. (O) « (Oq) + (O2) 
up to second order. Explicit expressions for jc and ps 
in terms of expectation values of boson fields are given 
in the Appendices. For clarity, we further divide the op- 



(0)\ 



erator O2 into two parts (O2) = {O. 
(O2*''') is calculated on the basis of L 



- (O^^^), where 
00 correlation 



functions, while {O^'') contains all finite length correc- 
tions. Then, analytical expressions for the charge cur- 
rent and spin density are readily obtained and given 



by (i 



2,c/ 



sgn {Vc, 



-Sgn (14) 



90- 



aa+- 



and 



respec- 



tively, where ujl — and Tcu is the short time cut- 
off. We note that the spin conserved scattering tc can 
only couple to Vc and, hence, reduce the total charge 
current in the voltage configuration of Fig. [D^a), while 
the spin-flip scattering t^, coupling to Vs, decrease the 
spin density in the configuration of Fig. [2][b). In the 
absence of finite length corrections, both charge current 
and spin density depend in a simple power law fashion 
on the applied voltage, in agreement with earlier work 
based on a renormalization group analysisi^ii^ii^ii^. Fur- 
thermore, we find that the ratio between charge current 
and spin density in the two dual voltage configurations 

is (ji^l) / (p'i^l) = s'^'/F' which can be used to obtain 
information about the scattering strength for the differ- 



ent types of scattering. In the finite length case, we use 
the numerical method of Ref. [l^l to evaluate (j2,c) and 
{p2,s)- As shown in Fig.lSja) and (b), the finite length wfll 
introduce oscillations in both the backscattering current 
and the spin density, which originates from Fabry-Perot- 
type interferences of plasmonic excitation^"^ . However, 
since the power law go + — 1 is always larger than 1 

for positive 50 1 (j^^c) ((/52°s)) ^i^^ increase rapidly with 
Vc {Vs) and dominate the oscillatory corrections for large 
Vc {Vs). We conclude that finite length corrections are 
not very important in this setup. 




eV,/fieoL 




FIG. 3: (Color online) (a) The backscattering current as a 
function of charge bias Vc generated by the spin conserved 
tc-term for the voltage configuration of Fig. [U^a). The blue 
dashed line corresponds to the charge current {j2'l), while 
the red solid line additionally includes the finite length cor- 

rections. The current unit is 



■ . (b) is similar to (a), 

except that the spin density (p'^l) and (p2,s) as a function of 
the spin bias Va are generated by the spin-flip ts-term for the 
voltage conflguration of Fig. [U^b). The spin density unit is 

2(2^9+1 

ft^ujLj, • (''^ Possible two-particle backscattering terms: spin 
conserved backscattering gi± and spin-flip backscattering gsf. 
(d) The voltage dependence of the correction to the spin den- 
sity {p2.s) due to the two-particle spin-flip term g^f. The spin 



density unit is ^(^^^,2^ 



In all expressions above, we 



iLL 



go — 0.7 and gs = 1. 



Long junction case.- Now we consider the opposite 
limit d ^ A/ as shown in Fig. [Ijc), which could be 
achieved by gradually narrowing the QSH sample into 
a ID wire experimentally. In the following, we an- 
alyze finite size effects related to d and assume that 
L ^ d, hence L — >■ c» is a reasonable approxima- 
tion. For simplicity, we model the LJ as a sTLL 



described by Ht 



Hti 



Ht2 with Hti 



H 



to 



4 



H!o and iJfo 



hv 

~2 



^ J dx [dx 



+ 



Hf, 



^^dx\{dJs? + jT{dxVs?\ for \x\ < d/2. Com- 

pared to Eq. ([3]), we find tfiat tlie cliarge sector re- 
mains uncfianged at tlie step from hTLL to sTLL, but 
the spin sector shows a stepwise variation of the in- 
teraction parameter as well as the velocity of spin ex- 
citations. Note that for finite spin-orbit coupling gs 
can be renormalized away from its non-interacting value 



1, see 



19.20 



Since the Hamiltonian in the charge 
sectors remains unchanged, the interaction between two 
hTLL (as described in H^q) will not affect the charge 
current (jc) at all. For the spin sector, besides H^q 



H2 



/J'rf dx [gii_ cos ( 
gij^-term 



87r(ps) + (7s/ cos (a/Stt^s)] where the 
related to spin conserved backscattering 
h.c. while the ggf-ievm is related 

to spin-flip backscattering V'l ^V'h ^V'L.tV'-R.t + h.c.^ see 
Fig. EJc). These are the most important perturbations 
to the spin sector in the absence of impurity scattering. 
Since the gij^-term conserves spin, it will not influence 
the spin density. Thus, we focus on the gsj-tcrm be- 
low. Up to second order perturbation theory, we ob- 
tain the following correction to the spin density (p2,s) — 

- %L%d lo /o^ '^^fi^^ ^) ^^•''"•3 9s = gs9o and 
the function f{r,R) is specified in Eq. (jB27|) of App. B. 
The above integration can be easily evaluated numeri- 
cally (similar to the Coulomb drag problem in^i) and 
the obtained correction to the spin density is shown in 
Fig. intd). To make analytical progress, we can again di- 
vide the obtained spin density into two parts (p2.s) — 

(P'il) + (/°2^i)' '^it^ (P-z^l) f^o^' ^^'^ infinite d and (P2^]) for 
the finite length corrections. It is found that (p'^l) al- 
ways dominates over (p'^l). In the limit eVs/hud ^ 1, 

-2 



we obtain {p. 



sgn {eVs 



_i9sjdrjilas 



r 



eVs 



(0)x 
2,s/ 

with LJd — J!J-^'^^ 

" gad 

Conclusions.- We have analyzed the charge current 
and spin density in a four-terminal setup based on two 
hTLL coupled to non-interacting electron reservoirs. Dif- 
ferent types of scattering mechanisms between the edges 
are taken into account, particularly short junctions and 
long junctions. It is shown that different power law 
dependencies as a function of bias voltages applied to 
the four terminals can be used to distinguish the scat- 
tering mechanisms. A simple duality relation between 
charge current and spin density has been discovered. 
Remarkably, all spin-related observables can be mea- 
sured by straightforward charge measurements in the 
four-terminal configuration. It is interesting to ask the 
question whether the spin density in such a setup can 
also be measured by other means (e.g. as a test of the 
model). Taking typical values for hvf = 3eV-A and 
eVs = 5meV, we find that the zero order spin density 
is about (po.s) ~ 5.3/im~^. This may be detected by 



state-of-the-art local Faraday/Kerr rotation^S. 
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Appendix A: Non-equilibrium Keldysh formalism of 
the spin-charge transport 

In this Appendix, we want to establish the general for- 
malism used for the calculation of spin-charge transport 
properties of our setup. The generating function of the 
system which is the starting point for all our calcula- 
tions is most conveniently written down in the helical 
basis. However, to make contact with the notation used 
in large parts of the main text, its representation in the 
spin-charge basis will be presented later on using Eq. (1) 
of the main text to express the non-chiral boson fields 
in terms of spin-charge fields. For our non-equilibrium 
transport calculations we use the Keldysh formalismi^Ji, 
in which all the field operators have a time variable re- 
siding on the Keldysh contour. This contour consists of 
the upper (forward) and the lower (backward) branch. 
In real time representation, this time dependence enters 
the formalism via a so called Keldysh index r/, where ip^ 
and 9^ denote the boson and the dual boson field on the 
upper branch (77 = +) and the lower branch (77 = — ), 
respectively. In the helical basis, the generating function 
is then given by 




a=± 



exp z / dtY^ [-^ imivl.ei] + H^wi,ei] 



+ Y.{j dxJ^{v)^l{v) + J dxJi{v)ei{v) 



(Al) 



where Sq is the action of the free Boson field at the edge 
a — ±, which reads^^ 



=Y.v dt dx -dt^Wl 



2g{x) 



(A2) 



Here we have set fi = 1. Hy which models the coupling 
to the leads can be written asi^ 



(A3) 
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where /i^ and are the generahzed chemical poten- 
tials coupling the helical liquids with the leads. Usually 
/Lt^ 7^ and /i^ = for an electric chemical potential. 
However, to keep the whole formalism as general as pos- 
sible, we consider both contributions. Ht describes inter 
edge coupling terms giving rise to different scattering pro- 
cesses. The concrete form of Ht depends on the detailed 
mechanism of the scattering process under investigation 
and will be discussed below. Throughout this work we 
treat as a perturbation and present the general for- 
malism up to the second order. The source terms J^v'a 
and J^6'2 introduced in Z[J] allow for the calculation of 
arbitrary correlation functions of interest by evaluating 
suitable functional derivatives of the generating function 
with respect to J at J = 0. 

Of course all calculations can equivalently be per- 
formed in the helical basis as well as in the spin-charge 
basis. However, we first transform the full Hamiltonian 
into the spin-charge basis since most of the quantities we 
want to compute are naturally written in terms of spin 
and charge fields. In the spin-charge basis, the generating 
function can be written as 



notations 



A 



\ (r) / 

V2Jl{r) ] 



r') 



/I 
1 



Vo 



\ 


-1 

1 / 



Ca = 



/C++(r,r') C+- 

CA-'Cr^r') Cl- 

Q++(r,r') Q+- 

VS^+(r,r') Ql' 



(r,r' 
(r,r' 
-(r,r' 
■(r,r') 



J-++(r,r') 



J--+(r,r') 



•^A 

TV 



2?++(r,r') V 



-(r,r') 2?- 



(A7) 



(A8) 



(A9) 

■(r,r')\ 
-(r,r') 
-(r,r') 
-(r,r')/ 
(AlO) 



Z{J\ = 



H> 



exp < I 



c(s) '''c(s) 



with 



A=C,5 



J-oo 



where 



Cr (r,r') 



(^^(r)^^, (r'))o, VT{vX) = 



Qr'(r,r') 



{^lirK (r'))o, 



and 



, {el{r)^l (r'))o = '''(r',r))t. With the 

compact definitions above, it is easy to check that the 
generating function can be written as 



(A4) 



X—c,s 



=EA(-l*IC-i4.^-FUrQ*A) 



Nz 

dtHt[4>c,*s 



dt i dx 



2g{xY ^^^^'^ 



5^ = 



r,=± 

2g{xy 



dx 



(A5) 



and 



riy — 



where /i^ 
and /x^ = 



E 



dtd'j 



Now let us introduce the following 



(All) 



Here the superscript T represents the matrix transpose. 
Note that the matrix product in the above equations in- 
cludes an integration over space and time. In order to 
separate the free part of the generating function from 
the terms to be treated perturbatively, we apply the fol- 
lowing shift to the vector of boson fields 

$a = $a-Aa[Ja], AA[JA]=iCAg^JA, (A12) 



A + ,6 
^A 



and explicitly Aa[Ja 

Al^[A] = --j= I dr' {C^{r,r')d,„l{r' 



with 



+J■f(r,r')a,/x^(r'))+^ J dr'(Cf (r, r') + r/Cf (r, r')) J^r') 
+i J dr'(J-f (r,r') +77^if(r,r'))JA'(r') (A13) 



[d,^,l{xWi + d,ni{x)9l\, (A6) Af[J, 



/27r 



dr' {Qf {ry)^.^^t{r') 
+V^{r,r')dy,{r')) +i J dr'(Qf (r, r') + 7?Q^(r, r')) (r') 
dr'(Pf (r, r') + rjV^ir, r'))4{r'). (A14) 
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The generating function then factorizes as follows 

Z[J] = Za,c[Jc]ZQAJs\Zt[J] (A15) 

where 

^o,a[Ja] -e-^Jlc.j,^ (A16) 

Ca - QCaQ^ - 

Cf(r,r') ^f(r,r')\ 



C«(r,r') Cf(r,r') J-«(r,r') F^{v,v' 
Qf(r,r') X>f(r,r' 
VQf(r,r') Qf(r,r') I?f(r,r') Pf(r,r')y 



(A17) 

with Cf = 0(t - tOdV'AW, ^A(r')])o, (r, r') = - 
t)([¥'A(r),^A(r')])o and Cf = {Wx{v),^x{v')})^, and 
similar definitions for the other correlation functions P, 
and Q. ^t['/] is given by 



Zt[J] = 



A -iV 4fC"^^ 

g 2 *A "-^A ^ 



Nz 



(A18) 



Next we need to relate the physical quantities to the 
generating function. The density is defined as (/5a) = 

^da:{ip\) and the current as (jx) = e^dt{ip\) = 

—evf\J^dx{0\)- As already mentioned, expectation val- 
ues like {(fix) and (Ox) can be conveniently calculated from 
the functional derivatives of the generating function Z[J] 

A direct calculation shows 



with respect to J'^ and J^. 
that 



i 1 5Z[J] 



SZo,x[Jx] 

Wr(r) +Z,[0]Wr(r) 



2Z[0] 5Jl{v) 
5Zt[J] 



j=o 



(A19) 



j=o 



with a similar expression holding for {9x)- From (jA19p . 
we find that (fx) can be decomposed into two parts: one 
is the zero order term coming from Zo,A[t/], and the other 
one is the scattering term coming from Zt[J]. Conse- 
quently, any physical quantity O can also be divided into 
two parts (O) = (Oo) + (O2). After a lengthy derivation, 
the density and the current can be expressed as 



(paW) = \ -dxi^x) = (po,A(r)) + (p2,A(r)),(A20) 



(Po,A(r)) - dr'dx ( Cf (r,r') J-f (r,r') ) 



(P2,A(r)) = i 



(A21) 



dr'a.Cf(r,r')(/r(r' 



I dr'9.J-f(r,r')(/A'(r'))^, (A22) 



OaW) = -e«M/-a,(0A) = Oo.A(r)) + 02.A(r)), 



(A23) 



Oo,A(r)) = ^1 dr'a, ( Q«(r,r') (r,r') 



(A24) 



{J2.x{r))^-^evf\l^ I dr'a.Qf(r,r')(/r(r'),^ 



levf 



^ J dr'd,V^{r,r'){fi{r'))^ (A25) 



where 

n 



rip 

IX 



<5v + 



and 



Nz ^ e -J 

means the average along the Keldysh contour with re- 
spect to the action J2x ^o,x — J dtHt[^c + Ac, + As]. 
Here we use <I> instead of $ to keep our notation simple. 
The expressions (jA20p to (jA25P are the central results 
for our approach to the calculation of charge current 
and spin density. 



Appendix B: Charge current and Spin density 

In this section, we calculate the charge current and 
spin density in our model. Now we consider the real- 
istic electric chemical potential (/x^ ^ and fi^ = 0) 
so that only /i^ and fi^ do not vanish. Let us choose 
the form of /i in each edge as dxfJ,± — —^J'l,±S{x + 
L/2) + fJ.2,±S{x — L/2), then we have dxHc = —f^i,cS{x + 
L/2) + fi2,cS{x - L/2) with Mi(2),c = Mi(2),+ + A*i(2) - 
and ^x^J'^s = ^^^i,sS{x + L/2) + ^2,sS{x — L/2) where 
/^i(2).s = Mi(2),+ ^ ^J'l(2).-■ Moreover we define the volt- 
age on the upper or lower edge as eVa = — l^2,a and 
correspondingly eVc = e (V+ + V-) /2 — {fii^c — A^2,c) /2 
and eV, =e{V+-V-) /2^ [tn^, - fi2^s) /2. ' 

From the expressions (|A21I) and (jA24p .the zeroth or- 
der of the charge current and spin density are easy to 
obtain, given the free correlation functions which will be 
discussed in some detail in the next section. The unper- 
turbed contributions yield 



2e2 



-Vc. 



{PO,s) = 



2e 
hv 



-Vs. 



(Bl) 



where we have recovered h. 



The influence of the scattering depends on its detailed 
mechanism. Here we discuss two opposite cases sepa- 
rately, the short junction and the long junction. For the 
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short junction, the scattering Hamiltonian is given by by 



Ht = J dx6{x — xq) 

+i,{i>l,^{x)^in{x) - i>l^{x)i>Lt{x) + h.c.)] .(B2) 



ifci^Dt = ^"^^ S{x' - xo) J dri5{xi - xo) 

m,cr 

^giV^S:<,(r')g-iv^H<,(ri)^ _^ g7rC^(ri,r')' 

(/c(^')).=0, 

(/r(^')>t = *^^Six' - xo) J dv,6{x, - xo) 



(B5) 
(B6) 



In the bosonization form, we have 



[- 



+ 



„7rC=(ri,r') 



(5(a;' — Xo) J dri6{xi — xo) 



H, 



dx5{x-xo) itcY V V sHV^^t) 2^me \ 



giy5F*a(r')g-iv^*<,(ri) 



+it,^r,''^r,''^sm{y^^a) 



(B3) 



(B7) 



{f!{x'))t = ^^<5(x' - xo) I dv,6{x, - xo) 



E 



mae 



imV^(A*(r')-A*(ri)) 



with tc(s) = ^^' ^'^ = V2{(pc + (J'Ps) and = V2(^s + 
a6s). Given this concrete form of Hf, the operators 
and read 



^giv^*a(r')g-i\Af*«(ri)^ ^ g7rC*(ri,r') 



(B8) 



= -it,V2^d{x' - xo)J2v'"'v 

COS V27T{(f'^ + ^^'^j') 

^ -iV2^<5(x' - Xo) 



,i<7 



X] V^^V^^o- COS V2^{<ft + 

a 

a<pt) + t,Y^ V'^V"^ cos (^+ + aOt) 

a 

ft{x') = -itsV2^S{x' - a;o) ^?7"^r?"^a 

a 

cosV2TT{ifj + a0f) 



where (• ■ •)t denotes the time ordered expectation 
value. In the latter equations we have used e^e^ = 

^A+B^i[A,B]^ ^g/^^ ^ gi(/^>o^ and (e*/(ri)e-'/(r2))o = 

(e-«/(ri)e«/(r2))p ^ e(/(ri)/(r.)-/^(0)>o. = ^(^^ + 

a^s), *a = V^((^« + aOs), C~ = (S^S^)o = 2 (Cc +C,) 
and C* = ($a$a)o = 2 [C, + + a (J", + Q,)], Ah,^ = 
^(^^ + (7^^) and A^,a = ^/2{Ai + aAl). 

As an example, we show the calculation for {fc{x'))t 
in some detail, 

(/c^(a;'))t « Y.'*c^S{x' - xo) f dnd{xi - xo) 

a ■' 

J2 V (cos (x/2^S+(r')) sin (\/2^S^(ri)) )^ (B9) 



^^^-^6{x' - Xo) J dri6{xi-xo) 



(B4) 



where ^1 = ^1+ A^'^^lO] and §1 = 91 + A^iO]. The 
average (A)* = (/^e-'^''''-/''^*^=[*+'^])o, is then given 



= E -^-^^i^' - Xo) J driS{xi - Xo) 

(7 



^g«V^=<'(r')e-»V^='<'(''i)^ + g7rCo(ri,r') 
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Then with the help of the expressions for the correla- 
tion functions listed in the next section, the second order 
correction to the spin density and charge current can be 
calculated as follows. 

(P2,s(a;)) =sgn(a;-2;o) / dti 



^^gimV^(A=(2;o,ti)-A=(3:o,0))g7rC=(xo,xo,ti) 



e 

and 



y^^g"nv^(^a(^o,*i)-^*(2;o,0))g7i-C*(2:o,a;o,ti) 

m.,a 
TrC^ {xo,xo,ti) 



Avf 



mae 



im^{Al{xo,ti)-Al{xo,0)) 



eti 



me 



(BIO) 



im^(A'^{xa,ti)~A'^{xo,'3)) 



(Bll) 



From the expressions for ^^(*^ hi the next section, we 
find that 



E 



(B12) 



and 



et^ r 



E 



^gimt 1 eVe g-n-C J (xo ,xo ,t 1 ) 



(B13) 



Numerical results obtained by the evaluation of the above 
integrals are presented in the main text. For the limit 
L — )■ (X), the correlation functions yield 



,1/1 

C^^{XQ,XQ,t) = -— go H 

27r \ go 



In 



C*(xo,xo,t) = [go + — 



In 



(B14) 



(B15) 



where lul = and Tcu is the short time cut-off. In this 
case, the correction to the spin density and the charge 
current can be calculated analytically to give 



90 + - 



P2,. 



TT^Tcu 



sgn(eT4) 



hwr 



(B16) 



and 



sgn (eVc) 



hjjT 



(B17) 



where h is recovered. 

Next we consider a long junction, which is modeled 
as a spinful Tomonaga Luttinger liquid. As shown in 
the main text, in the region x < \d/2\ where d is 
the length of the junction, the Hamiltonian is given 



H. 



to 



H, 



to 



Ht2 
2 



We find that H^^ = 
so that the Hamiltonian 



by Ht 

hv 
1 

of the charge sector has a constant interaction parameter 
go from — L/2 to +L/2. Therefore there is no scatter- 
ing in the charge sector. However, for the spin sector, not 
only the interaction parameter is changed across the junc- 

but there is 



tion (iJfo ^^Sdx {dj,f + ^{d^Vs) 



also an additional term H* 



19.20 



Hr 



dx 



gi± cos (^VSnip^^ + g^f cos (^VSnd. 



(B18) 



which comes from two particle scattering. In order to fo- 
cus on the new physics brought about by the long junc- 
tion during our perturbative calculation, we assume a 
hierarchy of length scales L 3> d ^ A/ which simplifies 
the problem. We can then take the L — ^ oo limit, which 
amounts to neglecting the influence of the leads. This is 
a reasonable assumption here because it is shown in the 
main text that already the finite size corrections with re- 
spect to d are small compared to the leading c? — oo 
contribution. Hence, finite size corrections with respect 
to L are negligible. Furthermore, we assume that no im- 
purities exist in the junction region, so that single particle 
backscattering can be neglected. This corresponds to the 
physical situation of a clean quantum wire. Within the 
above approximations, our model again becomes a g{x) 
model so that the formalism (jA20|) - (jA25|) can be applied. 
With the Hamiltonian (|B18|) . we have 

=gi±V8^sinV8^{ipt + At''') (B19) 
f! = g^fVS^ sin V8^{et + Ap<') (B20) 

Then up to the first order term we get 

(B21) 

.^g^V87e,(r')g-^V87e.(rl)^^ _^ g8,rI5, (n ,r' )1 (322) 



/gi\/87r¥>s(r')g-jV87r¥>3(ri)\ gSirCs (ri ,r') 
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It is convenient to divide p2,s into two parts 
(p2,s(r)) = (p2^.(r)) + (pL(r)) 



(B23) 



with 



(pLW) = |j§ J dx'dx,dt,^gn{x-x') 

m 

|^^gi\/87r<^a(£c',0)g-iV87ri^s(xi,ri)^_^ _ gSTrCs ,2;' ,ti ) 



2lo 



/ dx'(ia::idtisgn(a; — x') 
dd J 



E 



-j2m2^(xi-x')eya 

me "'f 



p8-7rCs(a:i,a:',ii) 



(B24) 



dx'dxidti 



j^g™y8¥(^^(^',0)-A^(a;i,ti)) 



dx'dxidti 



E 



(B25) 



Here we have taken /ii^s + /^2,s = and — gsSo- If 
we consider the region x > \d/2\, then it is easy to show 

that (jB24p vanishes because X^m "^s*^™ ^i)<=^s 
odd in r = — while Q^'^Co3{xifi,x ,t ) gygjj jj^ j.^ ]T;q_ 
(|B25|) can now be rewritten as 



{Pl.s) = 



■d/2 i-d/2 

Oik s\2^ I / dxif{x',xi) 

^{r^d) d J-d/2 J -d/2 



{gsfdfgs (B26) 
"'-l+l-RI 

— 22mri 4—1- 

me d 



4(fiw^)2d 

/■DO ^ 

87r-Do,(r,i?,ri) 



(B27) 



where r = , i? = ^^nd ri = iitj^ with 

uj^ — ji^. The fuh correlation function Vg {r,R,Ti) can 
be found in the next section. For comparison with the 
above single particle scattering process, we also divide 



P2,s into two parts P2,s = pIsI + P2,s2- The first part 
P2 si is the contribution for a quite long junction so that 
we can use the correlation function of the infinite wire, 

given by (r, R, r) = -j:^'^'^ (^ ^^"^^^ ^ "j which is 
independent of R. Then we can easily integrate over R. 
Furthermore all the singularities occurring during the in- 
tegrations appear in /92^g;^, so that our separate analysis of 
P2 si can help to keep the computing time of the whole 
calculation sustainable. In the infinite wire part P2sii 
/(r) is independent of R. Thus the expression (|B27|) can 
be simplified to 



{9sfdY~9. 



/I nOO T 7- 

dr dT'(-4)(l-|r|)sin(2j^T') 
-1 JQ l^d 



2(fia;^)2d J_ 
e^^^-^'i^y) sinSn^Vsir, r') 



(B28) 



This can be further simplified with the help of the fol- 
lowing equality 

dx {a + t{x + b))-" (a + i{x - b))-" e'P^ = 

O 

2_3/2 -ap /2|6|\ 

Here we require a > and fi > 1/2. With the above 
equality, Eq. (jB28p can be calculated to yield 



dr{l-\r\)\r\^-ilj^_,{\2^r\ 



eVs 



(B30) 



Similar to the calculation in Ref. [21], in the limit 
eVs/hujd 3> 1, the above expressions can be approximated 
as 



pLi ~ sgn{eVs) 



_2 4/5, 
/' / cu 



{gsfdf 



T^{2rgs){niol)U 



eVs 



(B31) 



Appendix C: Correlation functions 



In this section, we will provide some useful relations 
for the correlation functions used in the main text. First, 
we give a list of the zero temperature expressions for the 
correlation functions C, T ^ Q, and T) (derived in a similar 
way as in Ref. [T7|). 
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+i V In [""'u + + "^)']' ^ 1 (cu 



(C3) 



mez^,, xL.cu + (2? + m)2][r2„ + (2r? + m)2] 



(C4) 



/oo 
^^/giw(r-r')^ Q^(^ ?7 T — t') 
-oo 

short junction case, ^ = a;/L, r? = y/L, = ^ - ry, ^ , . 72e'2" sin(w^,) + sin(w^R) 

U = i + V,r = tujL, = ^ and r^u is the short = 2Sgn(Cr) (l - e^-lf-l) + i \_^-2g,L 

time cut-off. 5 and 7 = are different for charge and ^q^^ 
spin sector, g = = for the charge sector while „c» 

g = cjg = 1/go for the spin sector. For the long junction, / rfr'e*"^'^"'^ ^d^Q^i^^V^T ~ t') 

^ = x/d, r] = y/d, = ^ - V, = ^ + V, T = tuid and 

Tcu is the short time cut-off, where d is the length of the _ i ,j'^e'-^^' (•os(^'4,.) + c '~ vo^i^'^^u) 

junction. Now, = = and g = gc = 1 for the ~ 2^ ~'~ * 1 — ^2gj2ai 

charge sector, while = = ^ and g = g^ = g^g^ 
for the spin sector. 

Besides the above non-time ordered correlation func- 
tions, we have three other types of correlation functions: 
the retarded correlation function, the advanced correla- 
tion function, and the Keldysh correlation function de- 
fined in App. A. Useful identities for these correlation 
functions are 

/oo 
^^/gia,(r-r')5^C^($, r],T - /) = 

/ o^~2 i2cz, f-c ^^\ / dr' e'""^^-^'^ drJ''^{i,ri,T - t') 

^il fe^-IC.I + 2(7^6'^" cosjuj^r) + 7e'" cos{uj^r)) ^^ ' 



H l-^2,i2^ J _^...(, Aa:li.l\ I (7^e'^^sin((^g.)-V^sin(a;Cfl)) 

dr'e*^(---')a£C«(€, V,r-T')= ^ (^^^ 

.5" . _ ff(7^e'^-^sin(c^e.)+V'^sin(a;^fl)) j _Jr' e^''^^-^ ^ d^F^{^,n,r - r') 

(C6) 2 l-72e»2'^ ^ 
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/oo 
^^/giu)(r-r')g^piJ(^^ ?7, T - r') 
-CSO 

i 



whereas for spin part, only contributes, yielding 



/oo 
-OO 



(Cll) 



(C12) 

As seen in Sec. II, besides the correlation functions we 
also need A^'-^'(ri) — A^'-^'(r'), where A = c, s for our 
calculations, which are given by 

Anri)-^r(r') 

: J dV2 i^j^ dt3dt^C^{xi,tz,Y2) 



Al{v,)-Al{v') = - 



dT2 



(^J^^ dt3dt,Vf{xi,X2,t3-t2) 

+ J dX3dx^'D^{X3,X2,t' -t2)j dx^lJ.s{X2) 



2lT , 



(C17) 



dr2 



(rti rx, 
J dhdt,T^{x^,h,V2) + J 

as well as 

Ai{r,)-Ai{r') 



dxid^^F^{x3,t' ,V2) 



(C13) 



i^j^ dt:idt^Tf-{xi,X2,h - i2) 

(a«i,s - M2,s) - a;') . (C18) 



2\/27ra;£L 



= --j=jdr2 ^y dt3dt3Q^{xi,X2,t3 - 12) 

+ J dX3da;3Qx{X3,X2,t' -t2)j dx^lJ,l{X2) 

+ (^y «3p^(a;i,X2,i3-i2) 



(.C14) 



For the short junction, using the relation between S, ^ 
and ip, 9, we obtain 



Since for the charge part only n"^ contributes, we have 
At{r^)-At{r') 



H -X2\-\x' - X2\) I da;^llc{x2) 



1 



2V27r 
9. 



[{h - t')(Mi,c - A*2,c 



^E,a(a;o,ii) - A^^„{xo,t') = -^{h - t')eV^, (C19) 

A^,a{xo,ti) - A^^a{xo,t') = -^{h - t')eV,. (C20) 



(mi,c + M2,c)(a;i - x') 



(C15) 



A^(ri)-A^(r') 

dX2 



dx3 



2ojlLgc 



dx2l^c{X2) 



(Xi - x')(/ii,c - /U2,c 



2\/2'KgcU)1L 



For the long junction, we only need to change the param- 
(C16) eters L to d and lvl to uij,. 
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We assume the same parameters Vf and go at the two edges 

(a = ±). This is reasonable for two edges of the same 

system. 

We mention here the validity regime of our perturbation 
theory. Note that for go < 1 and gs ^ 1, the power law 
exponent — 2 will always be larger than 2. Thus, the 
second order correction to the spin density (p2,3) will in- 
crease faster than the zero order term (po,s) as a function 
of bias voltage. Hence, our perturbation theory is only jus- 

tified as long as re^/- (|^)' ||^|^"' « 1. 



